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Abstract 


The technique of computerized tomography (CT) has established itself as a leading 
tool in diagnostic radiobgy over the past thirty years and is catching on fast in the 
non-destructive evaluation area in a variety of situations. 

image reconstruction from projections are based either on transform methods 
or series expansion methods. Series expansion methods made their first appear- 
ance in the scientific literature and m tfre CT scanner industry around 1970. Series 
expansion methods are used when limited anrtount of data is available. Series 
expansion methods include Algebraic Reconstruction Technique and other opti- 
mization methods. Algebraic Reconstruction Techniques can broadly classified in 
two groups, Additive Algebraic Reconstruction Techniques (ART) and Multiplica- 
tive Algebraic Reconstruction Techniques (MART). 

Eight variations of ART and MART, are presented and compared in thb 
v«xk. Iterative tomographic algorithms have been applied to the reconstruction of 
simulated two-dimensional objects with internal defects. Comparison of the algo- 
rithms is based on their performance on the numerical generated objects. Impact 
of the variation of number of projections, projection angles, relaxation factor on 
the reconstructed images by these algorithms have been studied. 

The study shows that of the eight reconstruction algorithms applicable in the 
limited-data situation, those based on MART are fastest, most flexible and most 
accurate. The best of the algorithms is applied to the reconstruction of real data 
from 7 -ray scanner. This algorithm produces good reconstructions from limited 
amount of data. The data set includes projections at 4 angles over a 180° total 
angle of view (0°, 45°, 90°, 135°). 
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Nomenclature 


fj Exact value of the pixel 

fj Reconstructed field value of pixel 
Wij weight function 

(f)i exact projection of ray 

approximate projection of 
A relaxation factor 

e stopping criteria 

ai correction coefficient 

/^(pi difference in exact and approximate projection of ray 
C correction term 

k iteration number 

^max maximum value of weight 
N total number of pixels 

M total number of rays 

/ average of field values 

a average error 

d normalized RMS error 

r normalized absolute error 

E maximum error 

C normalizing factor 

N Number of photon counts after passing through the material 
Nq Initial number of photon counts 
/i(r, (p) Attenuation coefficient 



Chapter 1 


INTRODUCTION 


1.1 Reconstruction from projections 

Computerized Tomography (CT) is widely used in a large number of scientific, 
medical and technical fields. Hounsfield, inventor of first CT scanner, shared 
Nobel prize with A. M. Cormack in 1979, for their pioneering achievement in 
the field of computerized tomography. Algorithms for the image reconstruction 
from projection have extended our ability to visualize the internal structures 
of objects in a broad spectrum of physical applications ranging from molecular 
dimensions (electron microscopy) to cosmic dimensions (radio astronomy) . In 
words of Herman^^^ : 

’’Image Reconstruction from projections is the process of producing an 
image of two-dimensional distribution (usually of some physical property) from 
estimates of its line integrals along a finite number of lines of known locations” . 

In past, CT was done by X-rays but now other sources like 7-ray, acous- 
tic, laser, microwave, neutrons, seismic, magnetic resonance, radar etc. are 
also used. X-ray dose reduction, speed up of scan time, or other practical 
limitations which restrict the range of angles of projections or the density of 
rays in certain projections, all may lead to incomplete projection data collec- 
tion. Reconstruction from few directions arises in several nonmedical appli- 
cations. With limited data, iterative methods generally perform better than 
other techniques. The iterative method converts the reconstruction problems 
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into a system of linear eqxiations, with missing projectimi values treated like 
missing equations, i.e., they are ignored. 

In 7 -ray tomography, ray travels from source to receiver through the ob- 
ject. This ray contains the information of the part <rf object through which it 
passes. As the source and receiver are moved around ti» object, different pro- 
jections are obtained. By these projections, the distributitm ci some physical 
property like attenuation coefficients within the object is reomstructed using 
<me of the recoi^nKtion techniques. 


1.2 Reconstruction techniques 


Reconstruction techniques are broadly classified in two group; 

1. Transform methods or Direct methods. 

2. Series expansion methods. 

Direct methods include direct algebraic methods and Fourier Transform 
methods and Series expansion methods include Algebraic Reconstruction Tech- 
nique (ART), Multiplicative Algebraic Technique (MART), Simultaneous Iter- 
ative Reconstruction Technique (SIRT), Entropy Optimization etc. ART and 
MART are further classified into different classes depending upon the way 
corrections are applied to the algorithms. 

The basic difference between the series expansion methods and the trans- 
form method is that in the former case the field to be reconstructed is dis- 
cretized in the beginning, whereas in the later case continuous field is consid- 
ered and at end the final formulas are discretized for computational purpose. 

When the set of data available is large and complete, transform methods 
are generally used, whereas series expansion methods are preferred with in- 
complete data set. Series expansion methods are less efficient than transform 
methods (from computational speed and memory requirement point of view) 
but they can be used for irregular geometries. 

Using these reconstruction techniques, some physical property of object 
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is calculated like attenoiation coefficient or slowness (the inverse of acoustic ve- 
locity). In the ultrasonic tomography, the transit time &om source to receiver 
generally measxired and then tomographic image which is a distribution of 
slowness is reconstructed. This is also known as Ume-of-fli^ (TOP) tomog- 
raphy. 


1.3 Present work 


In the present work, diffierent series expansion methods, namdy ART and 
MART are studied. These reconstruction techniques are applied on the simu- 
lated data, for the reconstruction of a two-dimensional field. CosGauss phan- 
tom^ a numerically generated object, is used to compare various techniques 
in term of CPU time, convergence rates and errors. A situation in which the 
number of views is restricted but the total possible viewing angle is iK>t, i.e. 
projections at five angles equally spaced over a 180® total angle of view, is 
used for comparison of reconstruction al^rithms. The best algorithm among 
ART and the best among MART is then tested using a number of ima^ng 
geometries. The best among these two is further applied on various simulated 
two-dimensional fields which include Cosine phantom, circular plate with a 
central hole, square plate with five holes, square plate with two holes and 
square plate with five inclusions. This work has been carried out by Subbarao 
et but here, the following differences exist: 

• slightly different implementation of the algorithms, 

• different errors measrires are used, 

• fewer number of rays are used, consistent with the 7-ray CT scanner on 
which real data have been taken. 

Reconstruction is based on absorption coefficient data, as against simulated 
time-of-flight data used by Datta^^L 

The best reconstruction technique is then applied on real data, obtained by 7- 
ray CT scanner, to construct a two-dimensional field, which is the distribution 
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of abeorption coefl&den . Numerical computation were carried out cm a DEC 
400 workstation with 32 MB RAM and a 166 MHz processor. 


1.4 Thesis layout 

• Chapter 2 gives the detail of the reconstruction techniques, thear algo- 
liamsE various error measures used for aror anafysm of these aJgo- 
rithms- 

• Chapter 3 gives the detail of the experimental set-up used for sc annin g 
the specimen using the 7 -ray CT scanner. 

• Chapter 4 descnbes the application of the reconstruction techniques on 
the simulated fields. Different numerically generated fields used for test- 
ing the algorithms are discussed. 

• Ch^ter 5 provides discussions on the results based on the simulated and 
real data. 

• Chapter 6 summarizes the conclusions and suggests scope fear further 
works. 



Chapter 2 


THEORETICAL 

FORMULATION 


2.1 Formulation of a discretized model: 


A cartesian grid of square picture elements, called pixel, is introduced into the 
region that has to be reconstructed. These pixels are numbered following a 
regular fashion as shown in fig. 2.1. The attenuation coefficient is assumed to 
be constant say Xj throughout the pixel for j=l,2..,N. Source and detector 
are considered as points and the rays between them as lines. The length 
of intersection of ray with pixel, denoted by Wij for i=l,2...,M and 

j=l,2,...,N, represents the contribution of pixel to the total attenuation 
along the ray. 

The total attenuation of ray, denoted by <(>{, represents the line integral 
of the attenuation function along the path of ray. In this discretized model , 
the line integral takes the form of finite sum and the model may be described 
by a system of linear equations, as 

= (2.1) 

J=1 


In matrix notation 


<^=WX 


( 2 . 2 ) 
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DETECTOR 



Figure 2.1: Square pixels model for reconstruction from projection 

Problem now reduces to invert matrix W. The matrix W is quite sparse and 
its size may be large. Inversion of matrix W is quite impossible. Equation can 
be easily solved by series expansion methods. The goal of CT algorithms thus 
is to solve this set of M linear algebraic equations in N unknowns. 


2.2 CT using 7- rays 


Following Munshi^'^\ preliminary of 7-ray based CT is discussed below 
Single beam mono-energetic radiation attenuation phenomenon in a plane can 
be represented by 

I r I 

(2.3) 


N = NqBxp — J (f>)6l 


The value of fj. is characteristic of a material for a given type of radiation 
and the energy of radiation. Equ. 2.3 considers /i to be a two-dimensional 
function of position, as the path of radiation is assumed to be restricted to 
a plane. The energy dependence of p is also ignored. This is valid for 7- 
rays. Now, considering the material to have uniform attenuation coefficient. 
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equation 2.3 reduces to 

N = Noexp(-fiAl}, (2.4) 

In general, any cross-section of interest will have non-uniform distribution. 
Rewriting equation 2.3 as 


where, 


p = lni—) 


(2.5) 


( 2 . 6 ) 


Now, to recover p we require several p values measured along varioiis chord, c. 


2.3 Series Expansion Methods 


Different series expansion methods, namely ART and MART, are described in 
detail in the following subsections. 


2.3.1 ART 

Algebraic Reconstruction Technique (ART) was introduced by Gordon et 
for solving the problem of three dimensional reconstruction from projections 
in electron microscopy and radiology and was later recognized to be identical 
with Kaczmarz’s algorithm!-^^ for solving systems of linear equations. The first 
CT scanner developed by Hounsfield^^ used ART. Following Censor^^\ ART 
method used is presented below. 

An iterative process starts from an initial approximation f°eFCto the image 
vector. In an iterative step, the current iterate /* is corrected to a new iterate 
by taking into account only a single ray, say i*'*, and changing only the 
image values of the pixels which intersect this ray. The discrepancy between the 
measurement and the pseudo-projection data i=T,f=iWijfj) obtained 
from the current image /*' is redistributed among the pixels along the i*'* ray 
proportionally to their weights Wij in the whole ray. In this way, the pixel 
values along the ray are corrected to conform with the measurement 
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without changmg the rest of the image. Mathematically the whole process 
may be expressed as 


where, 


11^ Ir 

(2.7) 

JV 



(2.8) 

i=i 





(2.9) 

3=1 


||m*|p =< w\w* > 

(2.10) 


and 


Introducing the relaxation parameter, the modified step for updating the iter- 
ate can be expressed as 




( 2 . 11 ) 


where, A is the relaxation parameter, which is the real number, usually confined 
to the interval 


€i<A<2-e2, €i,€2>0. (2.12) 

If the correction is negative, it may happen that the calculated field of a recon- 
structed image becomes negative, in which case it is set to zero, as discussed 
by Gordon^^h 

ART may be further classified in different categories depending upon the 
correction applied. The four forms of ART are mentioned below. For aU these 
algorithms, let fj for j=l,2,..,N be the initial approximation of field value or 
pixel value. 

2.3.1. 1 Mayinger ART 

This is simplest iterative reconstruction algorithm proposed by Mayinger^^^h 
Pixel values are modified after all the correction values corresponding to the 
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individual rays are calculated. Let 4>i^ be the projection due to ray of view 
angle 6. Algorithm for Mayinger ART is presented below. 

For each view angle 9 


Step 1: For each ray ie, where ie denotes ray of view angle 9 and ie varies 
from 1 to Mfl, calculate 
(a) The approximate projection (0^,) as 


JV 


(b) Correction as 


4 = 

(2.13) 

i=i 


1 

II 

<1 

(2.14) 

function along each ray as 


N 



(2.15) 


i=i 


(d) Average value of correction as 

Mi, 




(2.16) 


Step 2: Repeat step 1 for all rays of view angle 9. 

Step 3: For each ray ie, modify the value of pixel, for all j=l,2,...,N, as 

/? = pi?) 

where A denotes the relaxation factor. 

Step 4: Repeat step 3 for all rays of view angle 9. 

Step 5: Repeat steps 1 to 4 for all view angles. This complete k*^ iteration. 
Step 6: Repeat steps 1 to 5 until 

fk-y-i _ fk 




X 100 < e 


(2.18) 


where e is the stopping criteria. 
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2.3. 1.2 Gordon ART 

In this methoS^^ correction is applied to all pixels through which ray passes, 
before correction for next ray is calculated. Thus, the number of rays in each 
view is not important. Algorithm is presented below. 

For each iteration k 


Step 1; For each ray i calculate 

(a) The approximate projection (^j) as 

(b) Correction as 

=4>i0- 4>it 

(c) Correction coefficient as 

i=i 


(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


(d) Modify the value of pixel, for all j=l,2,...,N, as 

XWijA<f>i 
Cti 


// = /r + 


( 2 . 22 ) 


where A denotes the relaxation factor. 

Step 2: Repeat steps (a) to (d) for all rays. This complete iteration. 
Step 3: Repeat steps 1 to 2 until 

X 100 < e (2.23) 

where e is the stopping criteria. 


2.3.1.3 Gilbert ART 

This algorithm^^^^ is also called as simultaneous iterative reconstruction tech- 
nique (SIRT). The pixel values are modified after all the correction values 


THEORETICAL FORMULATION 


11 


corresponding to individual rays are calculated. Thus, it is timilar to simple 
ART but correction applied to to the pixel value is different from it. Algorithm 
is presented below. 

For each iteration k 


Step 1: For each ray i calculate 

(a) The approximate projection (0,) as 

(b) Correction as 

— (pi (pi 

(c) Correction coefficient as 

ai = j:Wij^ 


(2.24) 

(2.25) 

(2.26) 


Step 2; Repeat steps (a) to (c) for all rays. 

Step 3: Identify all the rays passing through a pixel j and the corresponding 
values of Wij , A(pi and a,-. Let the total number of rays passing 
through pixel is denoted by Nj. 

Step 4: Modify the value of pixel, for all j=l,2,...,N, as 

(2.27) 

where A denotes the relaxation factor. This complete iteration. 
Step 5: Repeat steps 1 to 4 until 

X 100 < e (2.28) 

/*= 



where e is the stopping criteria. 
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2.3. 1.4 Anderson ART 

This algorithm^^^^ is also called as simultaneous algebraic reconstruction tech- 
nique (SART). SART combines both ART and SIRT. 

For each view angle 9 


Step 1: For each ray 1 $, calculate 

(a) The approximate projection (0^,) as 

i=i 

(b) Correction as 

(c) Correction coefficient as 

3=1 


(2.29) 


(2.30) 


(2.31) 


Step 2: Repeat steps (a) to (c) for all rays of view angle 9. 

Step 3: For each ray modify the value of pixel, for all j=l,2,...,N, as 


/> 

•'3 



+ 




(2.32) 


where A denotes the relaxation factor. 

Step 4: Repeat step 3 for all rays of view angle 9. 

Step 5: Repeat steps 1 to 4 for all view angles. This complete iteration. 
Step 6: Repeat steps 1 to 5 until 

ffc+i _ ffc 

i X 100 < e (2.33) 

fk - 


where e is the stopping criteria. 
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2.3.2 MART 

In ART, new field value is calculated by adding correction term to the initial 
value. While in multiplicative ART (MART), at each iteration, a new estimate 
of field value is calculated by multiplication of previous value by a correction 
term, i.e., 

= c X /;,i = 1, 2, ..., N. (2.34) 

where C is correction term, fj is the value of the pixel and k is the iter- 
ation number. Depending upon the correction term, MART may be further 
classified in four classes, as discussed by Verhoeven^^^^ Algorithm remains the 
same for all but correction term varies. Algorithm of MART may be stated as 

For each iteration k 

Step 1; For each ray i calculate 

(a) The approximate projection (0i) as 

(b) Modify the value of pixel, for all j=l,2,...,N, as 

Step 2; Repeat steps (a) and (b) for all rays. 

This complete iteration. 

Step 3: Repeat steps 1 and 2 until 

I pL X 100 < e (2.37) 

where e is the stopping criteria. 


(2.35) 

(2.36) 


Four classes of MART are 
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• GBH MART: This was proposed by Gordon et Correction term 
is given by 


C = 




(1 - ^) if 0 and Wij 0 


(2.38) 


1 else 

where A is relaxation factor and Wmax is the largest value of weight in 
W. 

• GH MART: This was proposed by Gordon and Herman^^^^. Correction 
term is given by 


1 — pp— ^(1 — if ^ 0 and Wij ^ 0 


1 else 


(2.39) 


where A is relaxation factor and Wmax is the largest value of weight in 
W. 

• Lent MART: This was proposed by Lent^^^K Correction term is given 

(2.40) 

( 1 else 

• Lent2 MART: A slight modification in Lent MART leads to Lent2 
MART^^^h Correction term is given by 

(^)AW-,-- if^.^0 
else 


-(? 


(2.41) 


W' ■ 

where Wij* is given by 


2.4 Convergence criteria 


It is necessary to determine when an iterative algorithm has converged to a 
solution which is optimal. In all algorithm, iterations are carried out until the 
following relation is satisfied 



100 <e 


(2.42) 


where, e is equal to 0.01. Here is the field value at {k + 1)*^ iteration and 
is the field value at iteration. 
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2.5 Error indicators 


Four error measures were used. They are used to find the deviation of recon- 
structed field from the exact field. Different algorithms may be compared on 
the basis of error measures. Lesser is the error induced by an algorithm in 
reconstructing a field, better is that algorithm. These can be applied only to 
the simulated field siace we know the exact field. CosGauss phantom is used 
for the error analysis of ART and MART. If fj is the exact value of pixel 
and fj is the reconstructed value of pixel, then the errors may be defined 
as 


1. AVERAGE ERROR: 


N 


2. NORMALIZED RMS ERROR: 


lE7=,[/y-/TJ 


(2.43) 


(2.44) 


where / is the average value of geometry used. The RMS error is large 
if there are few large errors in reconstruction. 


3. NORMALIZED ABSOLUTE ERROR: 

" EjLi |/;l 

This shows the effect of many small errors. 


(2.45) 


4. MAXIMUM ERROR: 

It is the largest difference between original and reconstructed field value. 
It may be defined as 


E = max I fj — fj I for j=l,2,...,N 


(2.46) 


Chapter 3 


Experimental setup 


The experimental unit consists of a collimator and a detector. Both are housed 
in the lead ring, whose inner diameter is 152 mm and outer diameter is 252 
mm. This unit can be used for collecting data for parallel beam geometry. The 
platform, on which specimen is placed, can have vertical movement of 105 mm 
and horizontal movement of 75 mm. It can also be rotated with a least count 
of 2°. Thus, this unit can have two translational and one rotational motions. 
The other equipments used are Nal scintillation detector, high voltage unit, 
low voltage unit, linear amplifier, single channel analyzer, timer, scalar. The 
source used is Experimental set up is shown in block diagram 3.1. 


3.1 Equipment 


The detail description of the equipment used in the experiment is given in 
the following sub-sections : 

3.1.1 Scintillation detector 

There are certain materials, e.g. organic, inorganic and plastic, 'which emit 
flashes of light or scintillations when radiation passes through them. These 
scintillations can be used to detect any type of radiations but are mostly used 
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for 7-ray counting. The output of a scintillator is proportional to the incident 
energy and hence can be used to find the energy of a 7-ray. 

Radiation is absorbed in the scintillator when it strikes the material. 
The absorbed energy is utilized in lifting the electrons to the excited state 
and part of it appeairs as heat energy. The electrons in the excited state de- 
excite in a very short time of about 10“®, leading to emission of light. This 
de-excitation time varies from material to material. 

The intensity of the emitted light is directly proportional to the energy 
lost by the incident radiation. The emitted light flashes are amplified by a 
photo multiplier. The photo multipilier tube consist of a photo cathode and 
several d5modes, each of which is maintained at a successively higher voltage. 
The photo cathode is usually a thin coating of Antimony-Cesium on the inside 
surface of the photo multiplier, whereas the dynodes are Antimony-Cesium, 
Silver-Magnesium or Copi>er-Beryllium. The electrons emitted by the photo 
cathode are amplified by releasing more electrons at each of the d5modes. The 
electrons from the last dynode are focused on to the anode connected to the 
preamplifier. The coupling of the scintillator to the photo multipHer is done 
by a light pipe. Normally, the scintillator is covered with a highly reflecting 
material, such as MgO and A/2O3, to avoid loss of light pulses. The light pipe 
is commonly made of glass, lucite or plexiglass. Its main purpose is to transmit 
efficiently the light pulses to photo cathode of photo multiplier tube. 

The advantages of scintillator for detecting 7-ray are: 

• Since all the light pulses are emitted in a very short time, the scintil- 
lator integrates these into a single pulse. The height of this pulse is 
proportional to the incident 7-ray energy. 

• Since the sdntiUator is solid, it has more atoms per unit volume than a 
gas-filled detector, leading to a better detection. 

• Scince the time constant of the scintillator is small, its dead-time is very 
small, facilitating a high rate. 

• A large scintillator increases the probability of absorbing the scattered 
7-rays, leading to an increase in photo peak efficiency. 
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Figure 3.1: Experimental set-up 
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A disadvantage in using a scintillator is the decrease in the light output with 
an increase in the energy of the incident 7-radiation. For example, the Nal 
(Tl) scintillator displays 20% variation in the light output per imit energy of 
7-ray over the range of 0-1 MeV. The light output also varies greatly with the 
typ)e of radiation. With the NeJ (Tl) scintillator, the electrons give a better 
response than alpha-particles of same energy. 

The overall efficiency of a scintillator is low because the detection pro- 
cess, namely, radiation causes light emmision which is then reconverted into 
electrons at the photo cathode, is complex. Typically 300 eV of energy is dis- 
sipated in the NeJ (Tl) scintillator for release of one photo electron, whereas 
in an organic scintillator about 1-5 KeV energy per electron is required. The 
corresponding energy for producing an ion-pair in a gas-filled datector is 30 eV 
and in a semiconductor is 4 eV. Scintillators are of many types, e.g. disk-type, 
right circular cylinder and well-typ>e. 

3.1.2 Linear amplifier 

It amplifies the output pulses from the detector. It has three gain controls, 
namely 

• Input attenuator, which have attenuation factors of 1,2,5,10,20 and 
50. 

• Coarse gain, which have amplification factors of 0.1, 0.3, 1. 

• Fine gain, with an amplification factors of 1-3. 

• Time constants, the pulse shaping provision of linear amplifier has 
variable integrating time constants of 0.1, 0.2, 0.5, 1, 2, 5 psec. 


3.1.3 Voltage units 

• High voltage unit: It supplies high voltage from the range of 300 to 
2500v. 
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• Low voltage unit : It supplies low voltage, ranging from -24 to 4-24. 

3.1.4 Single channel analyzer (SC A) 

An SCA can operate either in integrating (INT) or differentiating (DIFF) 
mode. In the INT mode, the lower level (base line) of the SCA acts as the 
discriminating bias. In the DIFF mode, the lower level is set at E and upper 
level at E' -I- AE, so that only the signal with a pulse height between E and 
E-fAE is counted by the scalar. The window AE can be continuously varied 
from 0 to 2 V. The base line adjustable between 0.2 to 10 V, thus making it 
possible to resolve and identify any pulse amplitude in this range. The output 
of SCA is constant and is of magnitude 10 V. 

3.1.5 Timer 

It sets the duration for which the counting is to be done. The factors are 1, 2, 
5 multiplied by 10 sec, 100 sec and 1000 sec. 

3.1.6 Scalar 

It is used to count the signals of pulse height between E and E+AE. 


3.2 Data collection 

Data collection using 7-ray CT scanner is done through following steps: 

• Identification of photo peak of the source 

Connect the low voltage and high voltage unit to the preamplifier of the 
scintillation head. Connect preamplifier to the linear amplifier. Connect 
the output from amplifier to single channel analyzer (SCA). Output from 
SCA is given to scalar through timer. Turn on the high voltage unit 
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and slowly increase its voltage to 750 V. Place the source through a 
collimator. Set window 'AE’ to 0.5 V and slowly increase the base from 
0.1 V. Set timer at 10 sec and collect the data at each base line setting. 
The counts will be higher at the photo peak of the source. 

• Collection of data 

Set the window and the base line at which photo peak was identified. 
Place the object between the source and detector on the platform. Set 
the timer at 100 sec and collect the data using parallel beam geometry 
collection mode. 

Setting of apparatus 

• Polarity^ -ve 

• Attenuator=20 

• Gain=l 

• Time constant=5 psec 

• Window width=0.b V 

• Base line=4.0 V 

• Mode= Differentiating 

• Timer = 100 sec 

• High voltage unit= 750 V 

3.3 Sources of error 

1. The detector must sense only those photons coming straight from the 
source. However, due to scattering some additional photons enter the 
detector, thereby increasing the count rate. Thus, density predicted by 
experimental results will be different from the true density of material. 
Radiation from these deflected photons is called secondary radiation. 
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2. Sources must be of small size in order to reduce scattering and geometric 
unsharpness. 

3. voltage fluctuation leads to shift in photo peak, or in other words there 
is fluctuation of Nq within a single pass. 

4. Lateral displacement of object by 1 mm, may not displace it exactly by 
1 mm due to gear mechanism. This may be true with rotational motion 
also. 

3.4 Precautions 

1. Counts must be taken for higher time period, to reduce error. 

2. To minimize error due to voltage fluctuation, each pass should be pre- 
ceded and succeeded by taking the value of Nq. If there is large variation 
in counts, the experiment for that case must be repeated else the average 
of the two values must be taken as the final value of A^o- 

3. To avoid backlash, object should be lateraly displaced and rotated in 
only one direction. 

4. The object should not be distrubed from its position till data collection 
completes. 

5. Source strength should be adequate for large No- 
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Validation of Programs 

4.1 Simulated data 


The programs for ART and MART have been tested on simulated data. Cos- 
Gauss phantom is used to compare different forms of ART and MART. 
CosGauss phantom is generated using the equation 


CosGauss{Xj y) 


where, 


C[0.3cos(x,j/) -I- 0 . 8 (e 2 :p{-[ 9 (ar - Xi)f - [ 6 ( 2 / - yi)]^} 
-f exp{-[{{x - a; 2 )]^ - [{y - ^ 2 )]^})] 

(4.1) 


cos(x, y) = 0.25{1 — cosf27r(a: + 

x{l — cos[27r(y -t- 0.5)^^^]} for j x |, | y |< 0.5 (4.2) 

= 0 else 


Thus, CosGauss phantom is the weighted superposition of the cosine 
phantom and two Gaussian peaks centered at xi=0.20, .yi=0.10, X 2 =- 0.20 
and y 2 =- 0 . 35 . Here, C is normalizing factor and is equal to 1.09. CosGauss 
phantom is shown in surface and contour plots in Figs. 4.1. The contour plot 
has the same contour levels. The contour interval is 0.1 units. 


A number of limited and unlimited data ima^ng geometries are used 
which are shown in Table 4.1, These may be basically divided into four cate- 
gories: 


1 . Limited angle with a limited number of views, in which both the 
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Figure 4.1: CosGauss phantom 


number of views and the possible viewing angle are restricted, e.g. cases 
A, B, C. 

2. Limited angle with unlimited number of views, where viewing 
angle is restricted but number of views are not, e.g. case D. 

3. Unlimited angle with a limited numbers of views, in which number 
of views are restricted but the total viewing angle is not, e.g. cases E, F. 

4. Unlimited angle with unlimited number of views, where both the 
number of views and possible viewing angle are not restricted, as in case 
I. 

The field to be reconstructed is discretized into 50x50 square grids of equal 
size. ART and MART performances are compared using case F, where, data 
sets consist of projections at five angles equally spaced over the viewing angle 
of 180'’ (0°, 36°, 72% 108°, 144°). In all 312 rays are considered for 5 projections. 
The ray location is chosen such that the ray spacing, perpendicular distance 
between rays, remains same for all views. Ray spacing equals to the dimension 
of a pixel. 
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Cases 

Projections 

No. of rays 

View angle 

A 

2 

100 

90 


3 

170 

90 

mm 

5 

294 

90 

D 

10 

616 

90 

E 

4 

240 

180 

F 

5 

312 

mKEHi 

G 

9 

566 

mHESlHi 

H 

10 

624 

180 

I 

^ 18 

1132 

180 


Table 4.1: Imaging geometries 


Table 4.2 compares the performance of ART family for relaxation factor 
0.5 while Table 4.3 compares the members of MART family at relaxation 
factor 1.0. GBH does not converge for higher values of relaxation factor. 
Oscillating pattern of GBH during iterative process is shown in Fig. 4.2 for 
relaxation parameter 0.01. Table 4.4 shows the performance of GBH at relax- 
ation parameter 0.001 for convergence criteria 0.05 and at relaxation parameter 
0.0001 for convergence criteria 0.01. 


5 Pro; 

ections, 312 rays, 180 view angle, A = 0.5 

Algorithnxs 

of 

ART 

Average 

error 

(%) 

Normalized 
RMS error 

(%) 

Normalized 
absolute 
error (%) 

Maximum 

error 

CPU 

time 

(0 

Iterations 

Mayinger 

10.54 

59.12 


0.4280 

8.13 

6 

Gordan 

10.11 

55.31 

53.31 

0.4084 

9.74 

8 

Gilbert 

6.90 

38.34 

36.40 

0.4113 

40.53 

39 

mssmsmi 

10.11 

55.32 

53.29 


9.31 

7 


Table 4.2: Performance of ART for CosGauss phantom 


The reconstructed images of CosGauss phantom by these eight algo- 
rithms are shown in Figs.4.3 to 4.10. 

Effect of relaxation parameter on errors is shown in Tables 4.5 and 4.6 
for Gilbert ART and Lent2, respectively. Here, relaxation parameter is var- 
ied from 0.01 to 1.99 and the corresponding variations in errors are studied. 
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5 Projections, 312 rays, 180 view angle, A = 1 

Algorithms 
of MART 

Average 

error 

(%) 

Normalized 
RMS error 

(%) 

Normalized 
absolute 
error (%) 

Maximum 

error 

CPU 

time 

(s) 

Iterations 

GBH 

- 


- 


“ 

- 

GH 

1.55 

11.00 

8.18 


1294 

1791 

Lent 

lllllllinni 

8.43 

6.69 

0.1083 

32L64 

432 

Lent2 

1.05 

7.51 

5.53 

0.1082 

114.37 

149 


Table 4.3: performance of MART for CosGauss phantom 


5 Projections, 312 rays, 180 view angle, A = 

0.001, conv 

ergence= 

0.05 

Algorithm 

Average 

error 

(%) 

Normalized 
RMS error 
(%) 

Normalized 
absolute 
error (%) 

Maximum 

error 

CPU 
time 
{ s) 

Iterations 

GBH 

3.73 

25.69 

19.65 

— 

0.3623 

— _ , 

241.36 

335 

jl 

11 

lections, 312 rays, 180 view angle, A = 

: 0.0001, con’ 

^ergence= 

=0.01 

Algorithm 

Average 

error 

(%) 

Normalized 
RMS error 
(%) 

Normalized 
absolute 
error (%) 

Maximum 

error 

CPU 

time 

(s) 

Iterations 

GBH 

3.55 

24.37 

18.73 

0.3214 

1293.6 

1780 


Table 4.4: performance of GBH for CossGauss phantom 



Figure 4.2: Oscillating pattern of GBH 
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CosGauss phantom is reconstructed vising Gilbert ART and Lent2 for 


Gilbert ART for 5 projections, view angle 180 

■ 


Average 
error (%) 

Normalized RMS 
error (%) 

Normalized 
absolute error (%) 

Maximum 

error 

0.01 

548 


46.12 

38.77 

0.5676 

0.1 

119 

6.62 

38.28 

34.90 

0.4583 

0.5 

39 

6.90 

38.34 

36.40 

0.4113 

1.0 

24 


43.22 

42.73 

0.4060 

1.5 

18 

12.88 

66.04 

67.93 

0.5042 

1.99 

3 

17.86 

87.69 

94.20 1 

0.5813 


Table 4.5: Variation of errors with relaxation parameter for Gilbert ART 


all the imaging geometries mentioned in Table 4.1. Tables 4.7 and 4.9 show 
the variation of CPU time and iterations with the number of projections for 
Gilbert ART and Lent2, respectively. Dependence of different errors, in the 
reconstructed field, on the number of projections is mentioned in Tables 4.8 
and 4.10 for these two algorithms. 

Reconstruction of CosGauss phantom, with projection sets of cases C 
and I, are shown in Figs. 4.11 and 4.12 for Gilbert ART and in Figs. 4.13 and 
4.14 for Lent2. 


lent2 for 5 projections, view angle 180 

A 

Iterations 

Average 
error (%) 

Normalized RMS 
error {%) 

Normalized 
absolute error (%) 

Maximum 

error 

0.01 

730 

1.29 

8.52 

6.78 

0.1062 

0.1 

179 

1.224 

8.14 

6.46 

■SlUifiH 

0.5 

250 

1.06 

7.32 


0.1104 

1.0 

149 

1.05 

7.51 

5.53 

■SUifiGH 

1.5 

213 

1.13 

8.54 1 

5.97 1 

0.1253 

1.99 

657 1 

1.06 

7.56 

5.60 

0.1104 


Table 4.6: Variation of errors with relaxation parameter for Lent2 


Lent2 is further applied to reconstruct physical objects with internal 
features. The field to be reconstructed is discretized in 50x50 pixels and 
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Gilbert ART for A = 0.5 

Case 

CPU time(s) 

Iterations 


6.53 

21 

B 

16.91 

32 

c 

32.64 

34 

D 

71.11 

35 

E 


36 

F 


39 

G 

76.06 

40 

H 

84.41 

40 

I 

155.89 

42 


Table 4.7: Variation of iterations and CPU time of Gilbert ART 


projection sets of cases C, F, I, are used for reconstruction. The different 
simulated objects used are: 


1. Circular plate with a central hole. This simulation is an object 
with absorption coefficient 1 everywhere outside the hole and 0 inside 
the hole. Circular plate is of diameter 50 units and holes have diameter 
of 20 units. Surface and contour plots of this simulated object are shown 
in Fig. 4.15. Reconstructed images are shown in Fig. 4.16. 


Gilbert ART for A = 0.5 

Case 

Average error 

{%) 

Normalized 
RMS error 

(%) 

Normalized 
absolute error 
(%) 

Maximum 

error 

A 

11.76 

69.79 

62.01 

0.5267 

B 

7.42 

45.21 

39.13 

0.4829 

C 

_ , 

42.96 

38.94 

0.4334 

D 

7.35 

42.56 

38.76 

0.4334 

E 

6.69 

36.76 

35.26 

0.3961 

F 

6.90 

38.34 

36.40 

0.4113 

G 

6.37 

1 34.10 

33.62 

0.3749 

H 

6.37 

34.08 

33.58 

0.3363 

I 

6.09 

31.92 

32.11 

0.2622 


Table 4.8: Errors in Gilbert ART 
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Lent2 for A = 1.0 

Cases 

CPU tmie(s) 

Iterations 

A 

3.96 

14 

B 

5.81 

12 

C 

139.55 

193 

D 

537.73 

366 

E 

9.34 

14 

F 

114.37 

149 

G 

238.81 

172 

H 

355.04 

233 

I 

1407.17 

522 


Table 4.9: Variation of Iterations and CPU time with projections in Lent2 


2. Square plate with five holes, this is simulated by considering absorp- 
tion coefficient 1 outside the holes and 0 inside the holes. Size of plate 
is 50 X 50 units and holes have diameter of 10 units. This is shown in 
surface and contour plots in Fig. 4.17. The reconstructed images are 
shown in Fig. 4.18. The plot of ll"‘, 26*^^ and 4r‘ rows of the exact and 
the reconstructed field for cases C, F, I are shown in Figs. 4.19 to 4.21. 

Here, the field values are either 0 or 1. The case in which the 
field values are widely separated is also studied. In which, holes have 
absorption coefficient 10 and plate has absorption coefficient 100. This 
case simulates a highly dense object with less denser materials embedded 
in it. The ability of algorithm to predict the presence of defects in denser 
media have been studied. The plots of ll"-. 26“ and 41" rows of exact 
field and the reconstructed field for cases C, F, I are shown in Figs. 4.22 

to 4.24. 

3. Square plate with two holes, this is similar to the object 2, except 
that it has only two holes. Object 2 is symmitrical while this object is 
non-symmetrical. This case have been considered to study the recon- 
struction of a non-symmetrical field. Fig. 4.25 shows the surface and 
contour plots of this object and the reconstructed fields m the form of 
contours src shown in Fig* 4,26. 
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Lent2 for A = 1.0 

Cases 

Average error 

(%) 

Normalized 
RMS error 

(%) 

Normalized 
absolute error 

(%) 

Maximum 

error 

A 

9.52 

65.62 

50.18 

0.4632 

B 

2.75 

17.90 

14.51 

0.1407 

c 

1.11 

8.02 

5.85 

0.1164 

D 

0.94 

6.89 

4.98 

0.0998 

E 

1.13 

7.20 

5.98 1 

0.0738 

F 

1.05 

7.51 

5.53 

0.1082 

G 

0.65 

4.64 

3.44 

0.0693 

H 

0.57 

4.17 

2.98 

0.0643 

I 

0.36 

2.57 

1.88 

0.0418 


Table 4.10: Errors in Lent2 


4. Square plate with five circular inclusions, in the above cases, the 
field values are either 0 or 1 which cleeirly distinguishes the holes from the 
object. This case have been used to see whether the algorithm is able to 
reconstruct a field , whose values are closer to each other. Here, square 
plate has absorption coefficient 1 and inclusions have slightly higher ab- 
sorption coefificient of 1.2. Fig. 4.27 shows the distribution of field values 
and the reconstructions are shown in Fig. 4.28. 

4.2 Experimental data 

The object, scanned by the 7-ray CT scanner, was made of plexi-glass, whose 
cross-sectional area was 2.5 cmx2.5 cm and was 1.5 cm thick. It has five holes, 
each of diameter 4 mm. Aim of the experiment is to reconstruct one plane of 
this object. It may be thus assumed as a square plate of size 2.5 cm by 2.5 cm 
with five holes. The field to be reconstructed is divided into 25x25 grids. The 
the ia 1 since the object is « horizontally 

by 1 mm. 

Object reas scanned by taking projections at 12 angles equaUy spaced 
over 180”. In ail 368 rays were taken. The object held (p) was reconstructed 

































Validation of Programs 


31 



Figure 4.3: Mayinger ART reconstruction of CosGauss phantom using case F 


using Lent2. The reconstructed field is the distribution of absorption coeffi- 
cient. 

Field is reconstructed for the following cases: 


1. Projections at two angles, i.e. 0°, 90°. Total rays are 48. 

2. Projections at four angles, i.e. 0°, 45°, 90°, 135°, and the total rays 
considered are 116. 

3. Projections at 12 angles and the total rays are 368. 

Reconstructed fields for the above three cases are shovm in Figs. 4.29 to 4.31. 
All the contour plots have the same contour levels, the contour interval is 
0.013 units. Plots for 6*'", 12*^^, 18*'^ rows of the reconstructed fields, by 4 
views and 12 views, are shown in Figs. 4.32 to 4.37. 


Validation of Programs 


32 



Validation of Programs 


33 



Figure 4.6: Anderson ART reconstruction of CosGauss phantom using case F 



Figure 4.7: GBH reconstruction of CosGauss phantom using case F 
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Figure 4.10: Lent2 reconstruction of CosGauss phantom using case F 



Figure 4.11: Gilbert ART reconstruction of CosGauss phantom using case C 
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Figure 4.17: Square plate with five holes 



(a) 



(b) 


(c) 


Figure 4.18: Reconstruction of square plate with five holes using (a) case C, 
(b) case F, (c) case I 
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Figure 4.19: Line plots of the 11^^ rows 


1.4 


1 \ 


1.2 



¥ 


“-r-jT ^ 

9 ^ / 

pi: / 

■jn'i ? 


exact — \- 
5(90) - 
18(180) -- 
5(180) - 


1 


1 0.8 
5 

2 
<D 

il 0.6 


0,4 


0.2 


0 L 
0 


d] i 




Field value pjQl^j yaiyg 


Validation of Programs 


41 



exact - 
5(90) - 
18(180) - 
5(180) -- 


Figure 4.21: Line plots of the 41** rows 
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Figure 4.22: Line plots of the 11** rows 
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exact — H 
5 ( 90 ) 
18 ( 180 ) 
5 ( 180 ) 


Figure 4.23: Line plots of the 26*'* rows 



exact - 
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18 ( 180 ) - 
5 ( 180 ) 


Figure 4.24: Line plots of the 41** rows 
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Figure 4.25: Square plate with two holes 






(b)- 


(c) 


Figure 4.26: Reconstruction of square plate with two holes using (a) case C 
(b) case F, (c) case I 
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5 
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Figure 4.29: Reconstruction with 2 views 




Figure 4.30: Reconstruction with 4 views 



Figure 4.31: Reconstruction with 12 views 


Experimental results 
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Chapter 5 


Results 


In this chapter , the performance of algorithms on simulated and real data is 
discussed. 


5.1 Simulated data 

5.1.1 Performance of ART 

Table 4.2 compares the performance of ART family. Mayinger ART requires 
fewer iterations to converge but errors in reconstructed field are higher among 
all the ARTs. Whereas, Gilbert ART takes higher CPU time to converge but 
the errors introduced in reconstructed field are lower. Gordon and .Anderson 
ARTs show nearly same performance. 

5.1.2 Performance of MART 

Table 4.3 compares the performance of MART family. GBH does no converge 
at higher value of relaxation parameter. Among the four M.ARTs, Lent2 is the 
best. It requires fewer iterations to converge and errors in the field predicted 
by Lent2 is lower. 
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GBH converges only for lower values of relaxation factor. It was tested for the 
relaxation factor of 0.01 and it produces oscillations during iterative process. 
It converged sharply to a value of about 0.7 within iterations about 40 and 
again diverged. It shows the convergence behaviour again after iterations of 
about 3000. GBH converged to 0.01 at very' low value of relaxation parame- 
ter of 0.0001 but requires a large number of iterations and the errors in the 
reconstructed field is maximum as compared to other MARTs. 


5.1.3 Effect of relaxation parameter on ART and MART 

Variation of errors and iterations w'ith relaxation factor is shown in Tables 4.5 
and 4.6 for Gilbert ART and Lent2 respectively. In Gilbert ART, as the value 
of relaxation factor is increased from 0.01 to 1.99, the iterations decrease but 
errors first decrease and then increase sharply'. 0.5 may be taken as the optimal 
value of relaxation factor for Gilbert ART. In Lent2, as the relaxation factor 
is increased from 0.01 to 1.99, the convergence time first reduces, become 
minimum at relaxation factor of 1.0 and then increases. Same is true with 
errors. Hence, for Lent2, best value of relaxation parameter is 1.0. 

5.1.4 Comparison between ART and MART 

Among the eight algorithms, MARTs are most accurate as compared to ARTs, 
although ARTs converge at lower iterations. Figures from 4.3 to 4.14 also ver- 
ify the better performance of MART over ART. Lent2 is superior among all, 
it reconstruct the field with lowest errors. The ability to determine 2500 un- 
knowns, with lower errors, from the set of 312 equations, shows the strength 
of Lent2. 

Gilbert ART and Lent2 is further tested on CosGauss phantom for the 
9 sets of projection data. Results are reported in Tables 4.7 to 4.10. 

• As the number of projections increase the error reduce, though conver- 
gence time increases. This is true, since more and more information is 



Results 


51 


now added or in simple words, more number of projections means more 
number of equations to solve 2500 unknowns. 

• Increasing the viewing angle decreases the errors. As tables 4.8 and 4.10 
predict that the field reconstructed by the data sets of viewing angle 90° 
is more erroneous than the field reconstructed by the data set of viewing 
angle 180°. 

• Orthogonality in the data sets also give the unique results. As seen from 
case E, which contains two sets of orthogonal projections at (0°, 90°) and 
(45°, 135°). Here data sets contains projections at 4 angles and total rays 
are only 240, but it shows better results as compared to cases C, D, F 
in the Gilbert ART, even though these cases contain more number of 
projections and rays. In Lent2, case E may not be considered superior 
to cases C, D, F but may be considered comparable with them. 

Figures also show that case 5(180°) is better than case 5(90°) and case 
18(180°) reconstruct the field which seems to be replica of the exact field. 

Lent2 begin the superior among all the algorithms is further applied to 
the 4 simulated objects. These objects are reconstructed using three important 
projection sets, that are cases C, F, I. 

The exact field and the reconstructed field are shown in the figures 4.15 
to 4.28. All the contour plots have the same contour contour levels. The 
contour interval is always 0.2 units. 

• Reconstruction by case C: Though it shows the location of holes or 
inclusions, but is not able to reconstruct them with exact shape and size. 
The circular holes take the form of ellipse. This is true, since the object 
is scanned from 0° to 90°. Plots of figures 4.19 to 4.21 also show the poor 
reconstruction of holes. The exact absorption coefl&cient of holes is 0 but 
the reconstructed value comes nearly 0.2, also the absorption coefficient 
of plate reduces to about 0.7. 

• Reconstruction by case F: Reconstruction improves as compared to 
case C. Still holes with exact shape and size are not reconstructed. Line 
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plots show the absorption coefficient of reconstructed holes and plate is 
about 0.2 and 1.2, respectively. Holes, although do not become perfectly 
circular but their shape improves as compared with case C. This also 
supports the argument that increasing the angle of view decreases the 
error levels. 

• Reconstructed by case I: This case gives surprisingly good results. 
Reconstructed field resembles the exact field. Line plots also show that 
plate and holes axe reconstructed with absorption coefficient 1 and 0 
respectively, except at the edges, where, a small error is present. These 
reconstructions also prove that larger are the number of projections, more 
accurate is the reconstruction. 

All the cases are able to predict the presence of inclusions in the plate. 
Although, reconstructions by cases C and F are not exact. Case I predict the 
exact location, shape and size of inclusions. This study shows the ability of 
algorithm to notice a small change in field values. 

Non-symmetrical case, i.e square with two holes, is also perfectly recon- 
structed. This result shows that the algorithm can be applied to any symmet- 
rical or non-symmetrical cases. 

If the field values are widely separated, i.e. the case in which absorption 
coefficient of holes is 10 and of plate is 100, algorithm is able to reconstruct 
the field. Here also, reconstruction by 5(180°) is better than 5(90°) and the 
reconstruction by 18(180°) is the best. Increasing the field values increases the 
error levels. 


5.2 Experimental data 


Field is reconstructed using 3 cases, which include reconstruction by 2 views, 
reconstruction by 4 views and reconstruction by 12 views. Figures 4.29 to 4.31 
show the reconstruction for all the cases. 
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• Reconstruction by 2 views: Reconstructed field is shown in contour 
plot of figure 4.29. This plot consists of a number of contours at same 
level. Although, it consists of contours at five locations of holes, but due 
to presence of other contours, this reconstruction may not be considered 
as perfect. 

• Reconstruction by 4 views: The reconstructed field consists of five 
holes, which can be clearly seen from the contour plot of figure 4.30. 
Plot shows the contours at the five hole locations. This result shows 
the uniqueness of orthogonal projections. The two sets of orthogonal 
projections , i.e (0°, 90°) and (45°, 135°), are able to reconstruct the 
field with very low errors. The ability of algorithm to determine 576 
unknown (pixel values) from the set of 116 equations (number of rays) 
using data sets which include projections at 4 angles (0°, 45°, 90°, 135°) 
is really surprising and the importance of this set of angles is evident. 

• Reconstruction by 12 views: Increasing the number of views must 
improve the reconstruction but here reconstructed field distorts. Holes 
become broad and take arbitrary shapes. It seems that, beside adding 
more informations, noise is indirectly added. Broadening of holes is due 
to shift in the information, i.e., reconstruction of a part by the informa- 
tion of some other part. This may be true if the lateral displacement of 
object by some distance, say x, does not move it exactly by that distance 
or due to improper rotation. Reconstructed field is shown in figure 4.31. 

Line plots of 6‘^, 12*^ and 18*^ rows of the reconstructed field by 4 views 
and 12 views are shown in figures 4.32 to 4.37. These figures, also indicate 
the broadening of holes. Broadening of holes in the field reconstructed by 12 
views is more than in the field reconstructed by 4 views. 



Chapter 6 


Conclusions 


In this study, four different MART algorithms and four different ART algo- 
rithms were programmed. These algorithms were compared with respect to 
their performance on CosGauss phantom, a numerically generated object. The 
following conclusions are drawn: 

1. Of the all algorithms tested, Lent2 gives best results in term of recon- 
struction accuracy. 

2. Among ART, Gilbert ART performs better, however its performance is 
poor when compared to MART. 

3. Relaxation parameter of 0.5 for Gilbert ART and 1.0 for Lent2 gives 
reasonably good results in terms of CPU time and errors. 

4. All the algorithms show a systematic behaviour with respect to the num- 
ber of projections and view angles. Increasing the number of projections 
decreases the errors in the reconstructed field but increases the CPU 
time. 

5. Increasing the angle of view decreases the error levels in the reconstructed 
field. 

6. Reconstruction improves when at least two projections are orthogonal to 
each other. 
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7. Inclusions present in the solid objects can be predicted reasonably well. 

8. Algorithms may be applied to any symmetrical or non-symmetrical ob- 
jects. 

Experimental results shows that, Lent2 produces superisingly good recon- 
struction from limited amount of data, even when the data has a significant 
amount of noise. The data set includes projections at 4 angles over a 180° 
total angle of view (0°, 45°, 90°, 135°). The uniqueness of this set of angles, 

i.e. two sets of orthogonal angles, is evident from this result. 

6.1 Suggestions for future work 

1. Dependence of the reconstructed field on the initial approximation of the 
field may be studied. 

2. The algorithms may be further tested on the numerically generated data 
with added noise. 

3. These algorithms may be applied to the ultrasonic TOF data, to test 
their performance on the non-linear data. 

4. 7 -ray CT scanner may be used to scan composite materials to study 
its ability to reconstruct the field of an object with varying absorption 
coefficient. 

5. Other series expansion methods may also be studied. 
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